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We study the texture of the exciton condensate at low temperatures in an independently gated 
electron-hole bilayer system. A model Hamiltonian is solved in real space within a mean-field approx- 
imation. It is found that, with increased electron-hole density polarization, the system experiences 
phase transformations from the zero center-of-mass momentum superfluid state, through one- and 
two-dimensional exciton pair modulated states, into the normal state. At weak density polarization, 
the modulating state resembles the Larkin-Ovchinikov state in superconductors in the presence of an 
exchange field in the weak-coupling BCS limit, and becomes stripe- like in the strong coupling BEC 
limit. In the one-dimensional modulated phase, the density of states exhibits low-energy intra-gap 
resonance quasiparticle states, which are localized in the nodal region. 
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The Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state 0, 
Q represents a class of unconventional superconducting 
states, where the superconducting order parameter mod- 
ulates in real space. It is dual to the other more fa- 
miliar class of unconventional superconducting states, 
in which the superconducting order parameter modu- 
lates (or even changes sign) in momentum space. In the 
FFLO state, the Cooper pairs form between two elec- 
trons with (k + q/2, f) and (— k + q/2, j.) of momentum 
and spin configuration. The center-of-mass momentum 
q is dependent on the extent to which the Fermi sur- 
face is Zeeman split by an exchange field. Although the 
FFLO state was predicted more than four decades ago, 
an undisputed verification of this state in superconduc- 
tors remains a great experimental challenge. A major 
reason lies in the fact that, in most superconductors, the 
Pauli paramagnetic effect of the applied magnetic field is 
negligibly small compared with the orbital breaking ef- 
fect. Owing to recent advances in the discovery of new 
materials and new technology, there is a revival of inter- 
est in the FFLO-type modulated states in a wide range of 
systems. Among them, newly discovered heavy fermion 
superconductors |3|,|4l, [5| , low-dimensional organic super- 
conductors Ja, hi ]8L 1 9. lOll. and trapped cold Fermion 
atoms p^l. [T2L I13I. I14II5I. Ilrj] are good candidates for the 
emergence of a superconducting FFLO state. In paral- 
lel, semiconductor bilayer systems [TtJ represent another 
promising context for the concept of FFLO-like states 
and the fundamental physics of the BCS-BEC crossover, 
in which there is condensation of electron and hole pairs. 

In electron-hole bilayers, the density of the electrons 
and holes can be varied independently of each other. 
The effect of electron-hole density imbalance here re- 
sembles the exchange effect of an applied magnetic field 
in a superconductor. Such an effect on the BCS-BEC 
crossover in semiconductor electron-hole bilayers has re- 
cently been investigated at zero temperature by Pieri and 
co-workers 18|, where an instability toward to the FFLO 



phase was identified. A more recent study |l9|. by treat- 
ing the FFLO phase and the Sarma phase [20[ on an equal 
footing, indicated that the FFLO state can be well sta- 
bilized by order parameter mixing effect. We note that, 
though collectively known as the FFLO state, the form 
proposed by Fulde and Ferrel (FF) [l| and that by Larkin 
and Ovchinnikov (LO) Q are slightly different. In the 
FF form, the order parameter has a homogeneous mag- 
nitude but with a modulated complex phase factor, that 
is a single- wavevector plane wave; while, in the LO form, 
the order parameter is real and spatially modulated. In 
Ref. 13, [l9|, the authors solve the problem directly in 



momentum space, which limits consideration to the FF 
form. In this paper, we consider a tight-binding model 
for an electron-hole bilayer and solve the BCS-type equa- 
tion self-consistcntly in real space. It is found that the 
LO form is always stabilized. Of particular interest, we 
find a crossover from the LO state to an exciton stripe 
state as the pairing interaction is tuned from the BCS 
limit into the the BEC limit. Our local density of states 
calculations show the existence of intra-gap states, which 
are localized on the nodal lines, in the one-dimensional 
LO phase. 

We start with a generic two-layer system, schemati- 
cally shown in Fig. Q] The two layers are separated by 
an insulator such that the interlayer tunneling is negligi- 
ble. The chemical potential in each layer can be tuned by 
gate voltages V 3i+ and V g -. In the tight-binding model 
for a square lattice in each layer, the Hamiltonian can be 
written as: 
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Here cj CT (c,:<j) are the creation (annihilation) operators of 
an electron at site i in layer a = 1 or 2. The quantities 
tij >(7 and fi a = Vg,cr are the hopping integrals and chem- 
ical potentials in each individual layer, respectively. We 
introduce gij as the inter-layer particle-particle coupling 
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FIG. 1: (Color) Schematic of a bilayer system. The carrier 
in each layer has electron or hole nature, which is determined 
by the individual band structure. The chemical potential in 
each layer is adjusted by the respective gate voltage. 

strength. In Hamiltonian [TJ the particle-particle inter- 
action within each individual layer is ignored, since they 
are expected to mainly renormalize each band. In our 
setting, explicit spin quantum numbers are omitted. As 
such, the index a can be regarded as a pseudospin index, 
which we will use intcrchangebly, a = +(— ) or f Q). 
In the mean-field approximation, the above Hamiltonian 
can be written as 

Hmf = ^[-%cr - ^kj]c\ a c ja + + h.c] 

ij,cr ij 

+ £|Ay| 2 /$«, (2) 

ij 

where the order parameter is defined as = fjij {cuc^} ■ 
Through a canonical transformation, the above Hamil- 
tonian can be diagonalizcd by solving the following 
Bogoliubov-de Gennes (BdG) equation: 

&)(?)-*(?)■ <3 > 

where the single-particle Hamiltonian and the pair po- 
tential are given by 

and 

A.4£<, r tanh(^) , (5) 

and (it™, w ™)Tr ans P° se are the eigenvector corresponding 
to the eigenenergy E n . The quasiparticle energy is mea- 
sured with respect to the Fermi energy. 

The above formalism is generic and applicable to both 
electron-electron or electron-hole bilayer systems, as de- 
termined by the band structure within each layer. Within 
the nearest-neighbor hopping approximation, the band 
dispersion for the clean case can be written as 

£ fe(T = -2t CT (cos k x + cos k y ) - ii„ , (6) 
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FIG. 2: (Color) Schematic of the energy dispersion for elec- 
tron (red) and hole (blue) layers. 

with t\ > and t% > for the electron-electron bilaycrs, 
and t\ > and < for the electron-hole bilayers. 
A schematic of normal-state band dispersion is shown 
in Fig. [21 Also the particle-particle interaction is long- 
ranged when the screening effect is weak, as considered in 
usual coupled quantum wells, and short-ranged when the 
screening effect is strong, which might be relevant when 
the metallic gate structures are present in the vicinity of 
quantum wells. For the purpose of the present discussion, 
we consider the electron-hole bilayer system, and consider 
a short-ranged Coulomb interaction with gij = gdij such 
that Aij = AiSij. 

We solve the above set of BdG equations self- 
consistently via exact diagonalization. Within the above 
mean-field treatment, the quasiparticle energy dispersion 
in the "quenched" uniform exciton condensate, that is, 
Aj = Aq, is of the form: 




FIG. 3: (Color) The spatial dependence of the exciton order 
parameter for changing h — (a), 0.2 (b), 0.6 (c), 0.7 (d), 1.0 
(e) and 1.60. The temperature is chosen as T = 10~ 3 . The 
average chemical potential is chosen as /i = —0.5. 
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From this energy dispersion form, one can see that the 
role of an exchange field in superconductors is played by 
the chemical potential /Ltx,2- However, the independent 
tunability of the hopping integral i 12 (or the effective 
bandwidth) within each layer has no counterpart in the 
case of superconductors. We note that the parameters 
ti t 2 are related to the effective mass mi^ appearing in 
the continuum limit 



IS Ha| . When £ k , 
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k,2: 



the 

first term on the RHS disappears, while the first term 
in the square is maximized. The quasiparticles are fully 
gapped, suggesting the "quenched" uniform state with 
zero momentum pairing is a truly favored state. How- 
ever, when £k.i 7^ — Ck.2, there emerges an instability 
toward a /mi£e-momentum pairing state. This analysis 
is consistent with that given in Refs. [HI, O and will 
be corroborated by the numerical calculations presented 
below. 

In our numerical calculations, we choose to measure all 
energy quantities in units of the hopping integral for the 
electron layer and set t\ = 1. The Boltzmann constant 
is set to ks = 1. A typical system size of Nl = 24 x 24 
is used throughout the calculation and the convergence 
criterion, set as the difference between two consecutive it- 
erations, is no larger than 3 x 10~ 5 . For the convenience 
of discussion, we introduce the average chemical potential 
and pscudospin polarization potential as fi = (/xi — 1^2)/^ 
and h = + /i2)/2, and the degree of bandwidth mis- 
match a = —fa/ti. Without loss of generality, we take 
a = I throughout this work. 

In Fig. [3l we show the order parameter texture on a 
square lattice for varying polarization field strength, h, 
at fixed values of fj, = —0.50 and go = 3. When the 
polarization degree of freedom h = 0, the exciton order 
parameter is uniform (see Fig. [3)Ja)). As h deviates from 
zero, the order parameter first exhibits a one-dimensional 
modulation (see Fig. |2Ib)). In this regime, the order 




FIG. 4: (Color) Schematic of the electron-hole bilayer phase 
diagram, covering three superfluid phases (USF, 1D-MP, and 
2D-MP), and several normal phases (ZZ, FF, ZF, FP, PP). 



parameter changes sign across the nodal lines. Our nu- 
merical calculations demonstrates (not shown) that the 
sign-change slope becomes steeper as the pairing strength 
is increased such that the nodal lines are replaced by 
sharp domain walls — a manifestation of a BCS-BEC 
crossover. With increased polarization, the modulation 
periodicity is shortened, and the one-dimensional mod- 
ulation gives way to a two-dimensional modulation (see 
Fig. EJc-e)). Simultaneously, the amplitude of the order 
parameter is decreased. Ultimately, the system enters a 
polarized normal state with even larger h (see Fig. 0(f)). 
Our calculations find that the LO-like state, which breaks 
the translational symmetry of the underlying lattice, is 
always favorable. 

In Fig. HI we provide a schematic of the typical \x-h 
phase diagram for the electron-hole bilayer defined on a 
square lattice. The phase diagram is constructed with 
delineating lines fi + h = ±4 and —fj, + h = ±4. For 
[i + h < —4, there are no electrons occupying the elec- 
tron band while for fj,+ h > 4, the electron band is fully 
occupied by electrons. Similarly, for —\i + h < —4, the 
hole band is fully occupied, while for —fi + h > 4 the hole 
band is empty. As such, we can map out the normal state 
phases, including a zero occupation phase (ZZ), full oc- 
cupation phase (FF), a phase with one of the layers fully 
occupied with the other empty (ZF), a fully polarized 
phase with partial occupation of only one type of carrier 
in the respective layer (FP), and a partially polarized 
phase (PP). The interesting exciton phases are limited 
in the area bound by the aforementioned four delineat- 
ing lines. They are categorized into a BCS (or BEC) 
uniform superfluid (USF) phase, and one- (ID) and two- 
dimensional (2D) modulated phases (MP). 

The inhomogeneous exciton state has physical impli- 
cations. We calculate the average local density of states 
for the electron and hole layers: 



Pe(E) 



1 J2\<\ 2 S(E~E n ), 



(8a) 



(8b) 



The Dirac 8 function in the above equations is approxi- 
mated by a Lorenzian 5(E — E n ) — > (j/n)[(E — E n ) 2 + 
7 2 ] -1 with 7 = 0.02. In the numerical calculation, we 
have also used a supercell technique to partially remove 
the finite size effect. In Fig. 03 wc show the average lo- 
cal density of states for various values of the polarization 
strength h. When h = (sec Fig. EJa)), the LDOS ex- 
hibits a standard BCS-like quasiparticle gap around zero 
energy, since the exciton condensate is uniform (corre- 
sponding to Fig. |3Ka)). In this phase, the quasiparticle 
density of states are identical for the electron and hole 
layers. When h is small but finite, for which the system 
is in the ID exciton modulated state (corresponding to 
Fig-d^b)), the density of states exhibits resonance peaks 
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FIG. 5: (Color) Average projected density of states for elec- 
tron (black) and hole (red) layers at various values of h = 
(a), 0.2 (b), and 0.7 (c). The average chemical potential 
H = -0.5. 





exhibited in the density of states for the ID modulated 
state, we present in Fig. [6] the spatial variation of the 
local density of states at the resonance energies. It can be 
seen that the local density of states at these energies has 
a maximum at the nodal lines and shows an exponential- 
like decaying behavior, indicating that these states are 
bound states. 

In conclusion, we have studied the texture of the ex- 
citon condensate at low temperatures in two parallel 
and independently gated electron-hole layers by solving a 
mean-field microscopic model Hamiltonian in real space. 
We have found that, with increased density imbalance, 
the system experiences transformations from the zero 
center-of-mass momentum supcrfluid state, through one- 
and two-dimensional exciton pair modulated states, into 
the normal state. At weak density polarization, the mod- 
ulated state resembles the Larkin-Ovchinikov state in su- 
perconductors in the presence of an exchange field in the 
weak-coupling BCS limit, and becomes stripe-like in the 
strong coupling BEC limit. In this phase, the density 
of states exhibits low-energy intra-gap resonance peaks. 
A local density of states imaging analysis indicates that 
these peaks correspond to bound state trapped near the 
nodal lines. Since the local density of states is propor- 
tional to the differential tunneling conductance as mea- 
sured by scanning tunneling spcctroscopy/microscopy, 
the inhomogencous exciton modulated state should be 
experimentally accessible 2l|. Therefore, an electron- 











hole bilayer system with density imbalance controlled by 
gate voltage, will provide an alternative setting for test- 
ing exotic phases proposed in superconductors. 
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FIG. 6: (Color) Spatial variation of the local density of states 
at E = 0.184 (a) and -0.190 (b), corresponding to the intra- 
gap peaks shown in Fig. 0L). 



within the BCS-like gap (see Fig.^b)). An equal number 
of such resonance peaks appear in the density of states 
for the electron and hole layers, respectively. We notice 
that in the BEC limit, multiple peaks can appear in the 
projected density of states. In addition, since the den- 
sity is unbalanced (for non-zero h), the overall density of 
states spectrum is shifted in opposite directions for the 
electron and hole layers. When h is further increased 
(see Fig. [5{c)), the shift of the overall density of states 
spectrum is increased for the electron and hole layers, 
and the quasiparticle gap-structure becomes invisible be- 
cause the exciton order parameter has been significantly 
suppressed (corresponding to Fig. 03d)). 

To understand the nature of the quasiparticle peaks 
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